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Abstract 

Imitative and contrarian behaviors are the two typical opposite attitudes of investors in stock markets. 
We introduce a simple model to investigate their interplay in a stock market where agents can take only 
two states, bullish or bearish. Each bullish (bearish) agent polls m "friends" and changes her opinion 
to bearish (bullish) if (1) at least mphb (mpbh) among the m agents inspected are bearish (bullish) 
or (2) at least mphh > mphb (mpbb > m-Pbh) among the m agents inspected are bullish (bearish). 
The condition (1) (resp. (2)) corresponds to imitative (resp. antagonistic) behavior. In the limit where 
the number TV of agents is infinite, the dynamics of the fraction of bullish agents is deterministic and 
exhibits chaotic behavior in a significant domain of the parameter space {phb, Pbh, Phh, Pbb,m}- A 
typical chaotic trajectory is characterized by intermittent phases of chaos, quasi-periodic behavior and 
super-exponentially growing bubbles followed by crashes. A typical bubble starts initially by growing 
at an exponential rate and then crosses over to a nonlinear power law growth rate leading to a finite- 
time singularity. The reinjection mechanism provided by the contrarian behavior introduces a finite-size 
effect, rounding off these singularities and leads to chaos. We document the main stylized facts of this 
model in the symmetric and asymmetric cases. This model is one of the rare agent-based models that give 
rise to interesting non-periodic complex dynamics in the "thermodynamic" limit (of an infinite number 
N of agents). We also discuss the case of a finite number of agents, which introduces an endogenous 
source of noise superimposed on the chaotic dynamics. 



"Human behavior is a main factor in how markets act. Indeed, sometimes markets act 
quickly, violently with little warning. [. . .] Ultimately, history tells us that there will be a 
correction of some significant dimension. I have no doubt that, human nature being what it 
is, that it is going to happen again and again." Alan Greenspan, Chairman of the Federal 
Reserve of the USA, before the Committee on Banking and Financial Services, U.S. House 
of Representatives, July 24, 1998. 

1 Introduction 

In recent economic and finance research, there is a growing interest in incorporating ideas from social 
sciences to account for the fact that markets reflect the thoughts, emotions, and actions of real people as 
opposed to the idealized economic investor whose behavior underlies the efficient market and random 
walk hypothesis. This was captured by the now famous pronouncement of Keynes (1936) that most in- 
vestors' decisions "can only be taken as a result of animal spirits - of a spontaneous urge to action rather 
than inaction, and not the outcome of a weighed average of benefits multiplied by the quantitative proba- 
bilities". A real investor may intend to be rational and may try to optimize his actions, but that rationality 
tends to be hampered by cognitive biases, emotional quirks, and social influences. "Behavioral finance" 
is a growing research field (Thaler (1993), De Bondt and Thaler (1995), Shefrin (2000), Shleifer (2000), 
Goldberg and von Nitzsch (2001)), which uses psychology, sociology, and other behavioral theories to 
attempt to explain the behavior of investors and money managers. The behavior of financial markets 
is thought to result from varying attitudes toward risk, the heterogeneity in the framing of information, 
from cognitive errors, self-control and lack thereof, from regret in financial decision-making, and from 
the influence of mass psychology. Assumptions about the frailty of human rationality and the acceptance 
of such drives as fear and greed are underlying the recipes developed over decades by so-called technical 
analysts. 

There is growing empirical evidence for the existence of herd or "crowd" behavior in speculative 
markets (Arthur (1987), Bikhchandani et al. (1992), Johansen et al. (1999, 2000), Orlean (1986, 1990, 
1992), Shiller (1984, 2000), Topol (1991), West (1988)). Herd behavior is often said to occur when 
many people take the same action, because some mimic the actions of others. Herding has been linked to 
many economic activities, such as investment recommendations (Graham and Dodd (1934), Scharfstein 
and Stein (1990)), price behavior of IPO's (Initial Public Offering) (Welch (1992)) fads and customs 
(Bikhchandani et al. (1992)), earnings forecasts (Trueman (1994)), corporate conservatism (Zwiebel 
(1995)) and delegated portfolio management (Maug and Naik (1995)). 

Here, we introduce arguably the simplest model capturing the interplay between mimetic and contrar- 
ian behavior in a population of N agents taking only two possible states, "bullish" or "bearish" (buying 
or selling). In the limit of an infinite number N —* oo of agents, the key variable which is the fraction 
p of bullish agents follows a chaotic deterministic dynamics on a subspace of positive measure in the 
parameter space. Before explaining and analyzing the model in subsequent sections, we compare it in 
three respects to standard theories of economic behavior. 

1. Since in the limit N — > oo, the model operates on a purely deterministic basis, it actually chal- 
lenges the purely external and unpredictable origin of market prices. Our model exploits the continuous 
mimicry of financial markets to show that the disordered and random aspect of the time series of prices 
can be in part explained not only by the advent of "random" news and events, but can also be generated 
by the behavior of the agents fixing the prices. 

In the limit N — ► oo, the dynamics of prices in our model is deterministic and derives from the 
theory of chaotic dynamical systems, which have the feature of exhibiting endogenously perturbed mo- 
tion. After the first papers on the theory of chaotic systems, such as Lorenz (1963), May (1976), (see, 
e.g., Collet-Eckmann (1980) for an early exposition), a series of economic papers dealt with models 
mostly of growth— Benhabib-Day (1981), Day (1982, 1983), Stutzer (1980). Later, a vast and varied 
number of fields of economics were analyzed in the light of the theory of chaos — Grandmont (1985, 
1987), Grandmont-Malgrange (1986). They extend from macro-economics — business cycles, models 
of class struggles, political economy — to micro-economics — models with overlapping generations, opti- 
mizing behavior — and touch subjects such as game theory and the theory of finance. The applicability of 
these theories has been thoroughly tested on the stock market prices — Brock et a/.(1987), Brock (1988), 
Brock-Dechert(1988), LeBaron (1988), Brock etal.{\99\), Hsieh (1989), Scheinkman-LeBaron (1989a, 
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1989b) — in studies which tried to detect signs of non-linear effects and to nail down the deterministic 
nature of these prices. While the theoretical models — Van Der Ploeg (1986), De Grauwe-Vansanten 
(1990), De Grauwe et al.(1993) — seem to agree on the relevance of chaotic deterministic dynamics, the 
empirical studies— Eckmann et a/.(1988), Hsieh-LeBaron (1988), Hsieh (1989, 1991, 1992), LeBaron 
(1988), Scheinkman-LeBaron (1989a,b) — are less clear-cut, mostly because of lack of sufficiently long 
time series (Eckmann-Ruelle (1992)), or, because the deterministic component of market behavior is 
necessarily overshadowed by the inevitable external effects. An additional source of "noise" is found to 
result from the finiteness of the number N of agents. For finite N, the deterministically chaotic dynamics 
of the price is replaced by a stochastic dynamics shadowing the corresponding trajectories obtained for 
N -> oo. 

The model presented here shows a mechanism of price fixing — decisions to buy or sell dictated 
by comparison with other agents — which is at the origin of an instability of prices. From one period 
to the next, and in the absence of information other than the anticipations of other agents, prices can 
continuously exhibit erratic behavior and never stabilize, without diverging. Thus, the model questions 
the fundamental hypothesis that equilibrium prices have to converge to the intrinsic value of an asset. 

2. We can also consider our model in the context of the increasing market volatility of financial 
markets. The volatility of prices generated by our chaotic model could give a beginning of an explanation 
of the excess volatility observed on financial markets — Grossman-Shiller (1981), Fama (1965), Flavin 
(1983), Shriller (1981), West (1988)— which traditional models, such as ARCH, try to incorporate (Engle 
(1982), Bollerslev et a/.(1991), Bollerslev (1987)). 

3. Finally, we can see speculative bubbles in our model as a natural consequence of mimetism. 
We can compare this to the two basic trends in explaining the problem of bubbles. The first makes 
reference to rational anticipations — Muth (1961) — and rests on the hypothesis of efficient markets. With 
fixed information, and knowing the dynamics of prices, the recurrence relation for the price is seen to 
depend on the fundamental value and a self-referential component, which tends to cause a deviation 
from the fundamental value: this is a speculative bubble — Blanchard- Watson (1982). This theory of 
rational speculative bubbles fails to explain the birth of such events, and even less their collapse, which it 
does not predict either. Recent developments improve on these traditional approached by combining the 
rational agents in the economy with irrational "noise" traders (Johansen et al. (1999, 2000), Sornette and 
A. Johansen (2001)). These noise traders are imitative investors who reside on an interaction network. 
Neighbors of an agent on this network can be viewed as the agent's friends or contacts, and an agent 
will incorporate his neighbors' views regarding the stock into his own view. These noise traders are 
responsible for triggering crashes. Sornette and Andersen (2001) develop a similar model in which the 
noise traders induce a nonlinear positive feedback in the stock price dynamics with an interplay between 
nonlinearity and multiplicative noise. The derived hyperbolic stochastic finite-time singularity formula 
transforms a Gaussian white noise into a rich time series possessing all the stylized facts of empirical 
prices, as well as accelerated speculative bubbles preceding crashes. 

The second trend purports to explain speculative bubbles by a limitation of rationality — Shiller ( 1 984, 
2000), West (1988), Topol (1991). It allows to incorporate notions which the neo-classical analysis does 
not take into account: asymmetry of information, inefficiency of prices, heterogeneity of anticipations — 
Grossman (1977), Grossman-Stiglitz (1980), Grossman (1981), Radner (1972, 1979). In our approach, 
which follows the second trend, the agents act without knowing the actual effect of their behavior: this 
contrasts the position of a model-builder — Orlean (1986, 1989, 1990, 1992). This, in turn, can lead to 
prices which disconnect from the fundamental indicators of economics. 

In the present paper we show that self -referred behavior in financial markets can generate chaos and 
speculative bubbles. They will be seen to be caused by mimetic behavior: bubbles will form due to 
imitative behavior and collapse when certain agents believe in the advent of a turn of trend, while they 
observe the behavior of their peers. 

Section 2 defines the model. Section 3 provides a qualitative understanding and analysis of its dy- 
namical properties. Section 4 extends it with a quantitative analysis of the phases of speculative bubbles 
in the symmetric case. Section 5 describes the statistical properties of the price returns derived from 
its dynamics in the symmetric case. Section 6 discusses the asymmetric case. Section 7 explores some 
effects introduced by the finiteness N < oo of the number of agents. In Section 8 we summarize our 
conclusions. 
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2 The model 



We consider an economy in which the population makes choices between two possible states when 
tomorrow's price is uncertain. The choice depends on expected capital gains. The portfolio choice 
then becomes a price expectation problem. Each agent has a different set of informations, obtained 
by observing other agents. Agents do not operate with reference to fundamental value, but rather with 
respect to expected market price. They are able to make profits if their expectations are judiciously 
chosen. It is rational for the agent — Keynes (1936), Orlean (1986, 1989), Sornette (2001) (see Chap. 
4) — to take into account collective judgments in order to make portfolio profits. 

That is why, for constituting expectations at time t + 1, the information used by an agent is the price 
expectation, at time t, of a certain sample of other agents randomly chosen among the population. This 
takes into account collective opinion and its expected correctness, that is, their confidence (or absence 
thereof) in the continuation of a deviation from the fundamental value. Their opinion refers to two kinds 
of price, market price and fundamental value, as exhibited by Keynes (1936): 

1 . Speculation relying on short term action and especially market opinion and market price. The 
most important aspect is the market price expectation, that is, the collective opinion about future market 
prices. 

2. Firm behavior: long term behavior relying on economic reality and fundamental value. This leads 
agents to detect excessive increase or decrease of market price and thus leads to anticipatory adaptation 
of the market price. This causes the collapse of the bubble. 

The importance of the interplay of these two classes of investing (which can be used by a same agent 
alternatively), corresponding to fundamental value investors and technical analysts (or trend followers), 
has been stressed by several recent works (Lux and Marchesi (1999), Farmer and Joshi, 2001) to be 
essential in order to retrieve the important stylized facts of stock market price statistics. This has recently 
been incorporated within a macroscopic model of the stock market with a competition between nonlinear 
trend-followers and nonlinear value investors (Ide and Sornette (2001), Sornette and Ide (2001)). We 
build on this insight and construct a very simple model of price dynamics, which puts emphasis on the 
fundamental nonlinear behavior of both classes of agents. 

These well-known principles generate different kinds of risks between which agents choose by arbi- 
trage. The former is a competing risk — Keynes (1936), Orlean (1989) — which leads agents to imitate the 
collective point of view since the market price includes it. Thus, it is assumed that Keynes' animal spirits 
may exist. More simply, there is the risk of mistaken expectation: agents believe in a price different from 
the market price. Keynes uses his famous beauty contest as a parable for stock markets. In order to 
predict the winner of beauty contest, objective beauty is not very important, but knowledge or prediction 
of others' prediction of beauty is. In Keynes' view, the optimal strategy is not to pick those faces the 
player thinks the prettiest, but those the other players are likely to think the average opinion will be, or 
those the other players will think the others will think the average opinion will be, or even further along 
this iterative loop. 

On the other hand, in the latter case, the emerging price is not necessarily in harmony with eco- 
nomic reality and fundamental value. Self -referred decisions and self-validation phenomena can then in- 
deed lead to speculative bubbles or sunspots — Azariadis (1981), Azariadis-Guesnerie (1982), Blanchard- 
Watson (1982), Jevons (1871), Kreps (1977). Thus, the latter risk is the result of precaution. It addresses 
the fitting of market price to fundamental value, and by extension, collapse of the speculative bubble. 

Both attitudes are likely to be important and are integrated in decision rules. Agents realize an 
arbitrage between the two kinds of risk we have described. That is why they have both a mimetic 
behavior and an antagonistic one: they either follow the collective point of view or they have reversed 
expectations. 

We are now going to put these assumptions into the simplest possible mathematical form. We assume 
that, at any given time t, the population is divided into two parts. Agents are explicitly differentiated as 
being bullish or bearish in proportion p t , and q t = 1 — p t , respectively. The first ones expect an increase 
of the price, while the bearish ones expect a decrease. The agents then form their opinion for time 
t + 1 by sampling the expectations of m other agents at time t, and modifying their own expectations 
accordingly. The number m of agents polled by a given agent to form her opinion at time t + 1 is the 
first important parameter in our model. 
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We then introduce threshold densities phb and phh- We assume < phb < Phh < 1- A bullish agent 
will change opinion if at least one of the following propositions is true: 

1. At least m ■ phb among the m agents inspected are bearish. 

2. At least m ■ phh among the m agents inspected are bullish. 

The first case corresponds to "following the crowd," while the second case corresponds to the "an- 
tagonistic behavior." The quantity phb is thus the threshold for a bullish agent ("haussier") to become 
bearish ("baissier") for mimetic reasons, and similarly, phh is the threshold for a bullish agent to become 
bearish because there are "too many" bullish agents. One reason for this behavior is, as we said, that the 
deviation of the market price from fundamental value is felt to be unsustainable. Another reason is that 
if many managers tell you that they are bullish, it is probable that they have large "long" positions in 
the market: they therefore tell you to buy, hoping to be able to unfold in part their position in favorable 
conditions with a good profit. 

The deviation of the threshold phb above the symmetric value 1/2 is a measure of the "stubbornness" 
(or "buy-and-hold" tendency) of the agent to keep her position. For phb = 1/2, the agent strictly 
endorses without delay the opinion of the majority and believes in any weak trend. This corresponds to a 
reversible dynamics. A value phb > 1/2 expresses a tendency towards conservatism: a large phb means 
that the agent will rarely change opinion. She is risk-adverse and would like to see an almost unanimity 
appearing before changing her mind. Her future behavior has thus a strong memory of her past position. 
Phb — 1/2 can be called the bullish "buy-and-hold" index. 

The deviation of the threshold phh below 1 quantifies the strength of disbelief of the agent in the 
sustainability of a speculative trend. For phh = 1, she always follows the crowd and is never contrarian. 
For phh close to 1/2, she has little faith in trend-following strategies and is closer to a fundamentalist, 
expecting the price to revert rapidly to its fundamental value. 1 — phh can be called the bullish reversal 
index. 

Putting the above rules into mathematical equations we see that the probability P for an agent who 
is bullish at time t to change his opinion at time t + 1 is: 

P = Prob {{x < m ■ (1 - phb)} U {x > m ■ p hh }) , (1) 

where x is the number of bullish agents found in the sample of m agents. 

In an entirely similar way, we introduce thresholds pbh, and pbb- The thresholds pbh and pbb have 
completely symmetric roles when the agent is initially bearish, pbh — 1/2 can be called the bearish 
"buy-and-hold" index. 1 — pbb can be called the bearish reversal index. The probability Q for a bearish 
agent at time t to become bullish at time t + 1 is: 

Q = Prob ({x < m ■ (1 - p bh )} U {x > m ■ p bb }) ■ 

We can combine these two rules into a dynamical law governing the time evolution of the popu- 
lations. Denoting p t the proportion of bullish agents in the population at time t, we can find the new 
proportion, pt+i, at time t + 1, by taking into account those agents which have changed opinion accord- 
ing to the deterministic law given above. To simplify notation, we let p t+ \ — p' and p t — p. Then, the 
above statements are easily used to express p' in terms of p, by using the probability of finding j bullish 
people among m (Corcos (1993)): 

P ' = p-p. Yl ( m )p m ~ j (i-py 

or j< m -{l-p hh ) 

+ (i-p). J2 ( m )(i-p) m_ v (2) 

or j«m (\-p bb ) 

= Fp. m {p) , 

where p = {phb, Pbh, Phh, Pbb}- Thus, the function F p . m (p) completely characterizes the dynamics of 
the proportion of bullish and bearish populations. 
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3 Qualitative analysis of the dynamical properties 



3.1 The limit m — ► oo 

The law given by Eq.^Jis not easy to analyze, and we give in Fig. [j] a few sample curves F pm . We see 
that as rn gets larger, the curves seem to tend to a limiting curve. Using this observation, our conceptual 
understanding of the dynamics can be drastically simplified if we consider the problem for a large number 
m of polled partners. Indeed, it is most convenient to first study the unrealistic problem m — oo and to 
view the large m case as a perturbation of this limiting case. The main ingredient in the study of the case 
m = oo is the Law of Large Numbers, which we use in a form given in Feller (1966): 
Lemma. Let g be a continuous function on [0, 1]. Then, for s £ [0, 1], 



3=0 V J 7 



ai(l-8) m -t-g(j/m) = g(s) . (3) 



We apply this lemma to the (piecewise continuous) function g = fh, where fh is the indicator 
function of the set defining P: 



_ ( 1 , if x > p M or x < 1 - phh, 
Jh \ , otherwise . 



(4) 



Similarly, we define 

f _ f 1 , if x > phh or x < 1 - p bb , 

Tb ~ \ , otherwise . { ' 

It is now easy to check that the lemma implies 

lim F p . m (p) = V -V Ml - P) + (1 - p) • fbip) = G p (p) . (6) 

m — >oo — — 

Note again that we do not consider G p (p) itself as an evolution law for the population of bullish agents, 
but G p can serve very well as an approximation for the true laws F p , m for large m. In Fig. [l] we show 
how the functions F p , m converge to G p . 



3.2 Classification of the different regimes 

In the preceding section, we have shown how to gain a qualitative understanding of the maps F PtTn , 
when m is large. We can now apply in a rather straightforward way the general theory of 1 -dimensional 
discrete time dynamical systems (see e.g., Collet-Eckmann (1980)) to the functions F p m . The recurrence 
Pt — > Pt+i can exhibit several typical behaviors which, for large m depend essentially only on the set 
of parameters p. We enumerate a few of them and refer the reader to Figs. || and||. In this section, we 
restrict our attention to the symmetric case phb — Phh and p^ = p^b- 

1. The most trivial case is the appearance of a stable fixed point. This will occur when the buy-and- 
hold index p^b — 1/2 is not too large and the reversal index 1 — phh is not too small. For example, this 
occurs for p^ = pbb = 0.75, phb = Pbh = 0.72, and rn = 60. Then, the population will equilibrate, 
and converge to p = 1 — q i=s 0.68, or to p ~ 0.32 (see upper panel of figure ||). 

2. The next more interesting case is the appearance of a limit cycle (of period 2): at successive times, 
the population of bullish and bearish agents will oscillate between two different values. This happens, 
e.g., for phh = Pbb = 0.76, with the other parameters as before (see second panel of figure |^). 

3. But for certain values of the parameters, e.g., p^h = Pbb = 0.85, the sequence of values of pt 
is a chaotic sequence, with positive Liapunov exponent (cf. Eckmann-Ruelle (1985)). The mechanism 
for this is really a combination of sufficiently strong buy-and-hold index phb — 1/2 and of sufficiently 
weak reversal index 1 — phh- This regime thus occurs when the opinion of a trader has a strong memory 
of her past positions and changes it only when a strong majority appears. This regime also requires 
a weak belief of the agent in fundamental valuation, as she will believe until very late that a strong 
bullish or bearish speculative trend is sustainable. Fundamentally, it is this self-referential behavior of 
the anticipations alone which is responsible for a deterministic, but seemingly erratic evolution of the 
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population of bullish and bearish agents. No external noise is needed to make this happen, and in general, 
we view external stimuli as acting on top of the intrinsic mechanism which we exhibit here (Eckmann 
(1981)). Note that the set of parameter values p for which chaos is expected (say, near the values used at 
the bottom of Fig. ^ has positive Lebesgue measure. 

We next consider in more detail the time evolution of p t for the parameter values of the last frame of 
Fig. & which are typical for the abundant set of "chaotic" parameter values, and we will show how the 
time evolution exhibits "speculative bubbles." This phenomenon is akin to the notion of intermittency 
(of "Type I") as known to physicists, see e.g., Manneville (1991) for an exposition. Indeed, we can 
distinguish two distinct behaviors in the last frame of Fig. |[ which occur repeatedly with more or 
less pronounced separation. The first process is the "laminar phase," which is seen to occur when the 
population p t is near 0.5. Then, the evolution of the population is slow, and the population grows 
slowly away from 0.5, either monotonically or through an oscillation of period 2, depending on p. This 
motion is slower when the inspected sample size (m) is larger, reflecting a more stable evolution for 
less independent agents. When the distance from 0.5 is large, erratic behavior sets in, which persists 
until the population reaches again a value of about 0.5, at which point the whole scenario repeats. The 
determinism of the model is reflected by "equal causes lead to equal effects," while its chaotic nature is 
reflected by the erratic length of the laminar periods, as well as of the bubbles of wild behavior. 

Having analyzed qualitatively the evolution of the number of bullish agents, we next describe how 
the price ir t +i of an asset at time t + 1 is related to the proportion p t of bullish agents. One can argue 
(Corcos (1993), Bouchaud and Cont (1998), Farmer (1998)) that the price change ir t +i — 7r t from one 
period to the next is a monotone function of p t (and, perhaps, of 7r t ). This function is positive when 
Pt > 1/2 and negative when p t < 1/2. If the reaction to a change in p t is reflected in the prices in the 
next period, then a bubble in p t will lead to a speculative bubble in the prices in the next period. Thus, 
our model predicts the occurrence of bubbles from the behavior of the agents alone. Furthermore, for 
quite general laws of the form 

7T t+ l = H(TT t ,p t ) , (7) 

a simple application of the chain rule of differentiation leads to the observation that the variable ir t has 
the same Liapunov exponent as p t . In fact, this will be the case if < d-^H < A and d p H > c > 0, 
where A is the Liapunov exponent for p t , as follows from <5"7r f+ i = d w H ■ Snt + d p H ■ 5pt. This condition 
is, in particular, satisfied for a law of the form Tit+i = tt* + G(pt), where G is strictly monotone. Thus, 
chaotic behavior of bullish agents leads to chaotic behavior of prices. 

In the sequel, we shall take the simplest form of a log-difference of the price linearly proportional to 
the order unbalance (Farmer (1998)), leading to 

lri7r t+ i - hnr t +i = r t +i = j(pt - 5) , (8) 

showing that the return r t calculated over one period is proportional to the imbalance p t — \. Thus, the 
properties of the return time series can be derived directly from those of p t as we document below. 

To summarize this qualitative analysis of the case of an infinite number N of agents, we observe a 
time evolution which, while satisfying certain criteria of randomness (such as possessing an absolutely 
continuous invariant measure and exhibiting a positive Liapunov exponent — cf. Eckmann-Ruelle (1985)) 
at the same time exhibits some regularities on short time scales, since it is deterministic. Our model thus 
establishes that straightforward fundamental conditions may suffice to generate chaotic stock market 
behavior, depending on the parameter values. If the market adjusts present market price on the basis 
of expectations and mimicry — self -referred behavior — then chaotic evolution of the population will also 
imply chaotic evolution of prices. 
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4 Quantitative analysis of the speculative bubbles within the chaotic 
regime in the symmetric case 

For an infinite number N of agents and in the symmetric case phb — Pbh = Pi and phh = Pbb = p2, let 
us rewrite the dynamical evolution (||) of the system as 

P'=P-pjt (7)^" (1 " pYf (£) +il - p) £ (7) (1 - Pr ' jpJf (m) ' (9) 



where 



Let us define 



which yields 



1 0, otherwise. v ' 



J=0 



p = F m (p) = p - p ■ g m (p) + (1 - p) ■ .9m(l - p) ■ (12) 

This expression (|l2|) generalizes (^J) to arbitrary m. 

As was described in the previous section, this system can exhibit chaotic behavior for certain values 
of the parameters. An example is given in figure Qwhich shows a long time series, showing many positive 
bubbles and negative bubbles interrupted by chaotic oscillatory phases. For the time being, we do not 
worry about the existence of the negative bubbles, which are rarely if ever observed in real markets: this 
is an artifact of the symmetry phb = pbh = Pi an d Phh = Pbb = P2, that we shall relax later. Keeping 
the symmetry assumption simplified the theoretical analysis without changing the key results obtained 
below. 

Let us consider the first bubble developing in the time interval from t — 35 to t = 546 as seen in 
figure ||-a). Figure ||-b plots the logarithm of p — 1/2 as a function of linear time: the linear trend from 
t = 35 to t w 480 seen in the lower panel qualifies an exponential growth p — 1/2 oc e Kt (with k > 0) 
followed by a super-exponential growth accelerating so much as to give the impression of reaching a 
singularity in finite-time. 

To understand this phenomenon, we plot the logarithm of F m (p) — p versus the logarithm of p — 1/2 
in figure^ for three different values of m = 30, 60 and 100. Two regimes can be observed. 

1. For small p — 1/2, the slope of \og 10 (F m (p) — p) versus log 10 (p — 1/2) is 1, i.e, 



p'-p = F m {p)~p~a{m)\p--j. (13) 

This expression ( |l3| ) explains the exponential growth observed at early times in figure ^| 

For larger p — 1/2, the slope of \og 1Q (F m (p) — p) versus log 10 (P — 1/2) increases above 1 and 
stabilizes to a value p(m) before decreasing again due to the reinjection produced by the contrarian 
mechanism. The interval in p — 1/2 in which the slope is approximately stabilized at the value 
p(m) enables us to write 



F m (p)-p~ P(m) [P-2J ' with V > 1 ■ (14) 
These two regimes can be summarized in the following phenomenological expression for F m (p): 

F m {p) = \ + (1 - 2 5ro (l/2) - <4(V2)) [p~ 5) +Km)(p-~y ( ] , (15) 
1 / 1\ / 1\ /l x 



. .p- 2 I +a(m) lp-- 1 +/9(m) (p-- J with/i>l, (16) 
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and 

a(m) = -2. 9m (l/2) - g' m (l/2) . (17) 

This expression can be obtained as an approximation of the exact expansion derived in the Appendix. 
In order to check the hypothesis (|l6), we numerically solve the following problem 



mm 



.(p)-5-[l + a] lP-- 2 )-P[P-~ 2 



(18) 



which amounts to constructing the best approximation of the exact map F m (p) in terms of an effective 
power law acceleration (see ( p0| ) below). The results obtained for m = 60 interacting agents and p^ = 
pbh = 0.72 and phh = Pbb = 0.85 are given in table [j] and shown in figure^. The numerical values of a 
are in good agreement with the theoretical prediction : a(m) = F^ n (l/2)—l which yields a(m) ~ 0.011 
in the present case (m = 60, phb = Pbh = 0.72 and phh = Pbb = 0.85). As a first approximation, we can 
consider that the exponent p is fixed over the interval of interest, which is reasonable according to the 
very good quality of the fits shown in figure ^. We can conclude from this numerical investigation that 
p(m) £ [3, 4]. A finer analysis shows however that the exponent p is in fact not perfectly constant but 
shifts slowly from about 3 to 4 as p increases. This should be expected as the function F m (p) contains 
many higher-order terms. We can also note that the parameter p c = {[3 /a) 1 ' M , which defines the 
typical scale of the crossover remains constant and equal to p c ~ 0.70 for all the fits (except for the 
largest interval p — 1/2 < 0.2, for which p c = 0.8). In sum, the procedure ( |l8l ) and its results show 
that the effective power law representation (|l6|) is a cross-over phenomenon: it is not dominated by the 
"critical" value phb = Pbh of the jump of the map obtained in the limit of large m. 

Introducing the notation e = p — 1/2, the dynamics associated with the effective map ( [T^ ) can be 
rewritten 

e' - e = a(m)e + P(m)e"( m \ (19) 
which, in the continuous time limit, yields 

= a(m)e + /3(m)e^ m) . (20) 

at 

Thus, for small e, we obtain an exponential growth rate 

e t ~ e a (™)* , (21) 

while for large enough e 

e t - (t c ~ t)~M^ . (22) 

For example, for m — 60 with phb = Pbh = 0.72 and phh = Pbb = 0.85, we can check on figure || 
that p(m) = 3, which yields for large e: 

pt -^vi=f (23) 



The prediction (23) implies that plotting (p t — l/2)~ 2 as a function of t should be a straight line in 
this regime. This non-parametric test is checked in figure M on five successive bubbles. This provides 
a confirmation of the effective power law representation ( |lq ) of the map. The fact that it is the lowest 



Optimization Domain 


a (3 /j, 


< p- | < 0.05 
<p - | < 0.10 
< p- | < 0.15 
<p - | < 0.20 


0.011 11.67 3.27 
0.013 43.66 3.77 
0.014 60.32 3.91 
0.004 30.64 3.54 



Table 1: Optimized parameters a, (3 and pi for several optimization interval with m=60 interacting agents 
and p hb = p bh = 0.72 and p hh = p bb = 0.85 
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estimate /j, w 3 shown in table yj which dominates in figure |] results simply from the fact that it is the 
longest transient corresponding to the regime where p is closest to the unstable fixed point 1/2. This is 
visualized in figure ||by the horizontal dashed lines indicating the levels p — 1/2 = 0.05, 0.01 and 0.2. 
This demonstrates that most of the visited values are close to the unstable fixed point, which determines 
the effective value of the nonlinear exponent p w 3. 

With the price dynamics (|j), the prediction ( p2| ) implies that the returns r t should increase in an 
accelerating super-exponential fashion at the end of a bubble, leading to a price trajectory 

7T f =7T c -C{t c -t)T^ , (24) 

where tt c is the culminating price of the bubble reached at t — t c when p(m) > 2, such the finite-time 
singularity in r t gives rise only to an infinite slope of the price trajectory. The behavior (^) with an 
exponent < ^[^l-i < 1 has been documented in many bubbles (Sornette et al. (1996), Johansen 
et al. (1999, 2000), Johansen and Sornette (1999, 2000), Sornette and Johansen (2001), Sornette and 
Andersen (2001), Sornette (2001)). The case to = 60 with phb = Pbh = 0.72 and phh = Pbb = 0.85 
shown in figure H leads to = 1/2, which is in reasonable agreement with previously reported 

values. 

Interpreted within the present model, the exponent of the price singularity gives an estimation 

of the "connectivity" number m through the dependence of /i on m documented in figure |6| Such a 
relationship has already been argued by Johansen et al., (2000) at a phenomenological level using a 
mean-field equation in which the exponent is directly related to the number of connections to a given 
agent. 



5 Statistical properties of price returns in the symmetric case 

Using the price dynamics (§), the distribution of p — 1/2 is the same as the distribution of returns, which 
is the first statistical property analyzed in econometric work (Campbell et al. (1997), Lo and MacKinlay 
(1999), Lux (1996), Pagan (1996), Plerou et al (1999), Laherrere and Sornette (1998)). Note that the 
distribution of p — 1/2 is nothing but the invariant measure of the chaotic map p'(p) which can be shown 
to be continuous with respect to the Lebesgue measure (Eckmann and Ruelle (1985)). Figure ^ shows 
the cumulative distribution of r t oc p t — 1/2. Notice the two breaks at \p — 1/2| = 0.28, which are due to 
the existence of weakly unstable periodic orbits corresponding to a transient oscillation between bullish 
and bearish states. 



Figure 10 plots in double logarithmic scales the survival distribution of r t oc p t — 1/2 for to = 30, 60 
and 100. For m = 60, we can observe an approximate power law tail but the exponent is smaller than 1 in 
contradiction with the empirical evidence which suggests a tail of the survival probability with exponents 
3 — 5. In the other cases, we cannot conclude on the existence of a power law regime, but it is obvious 
that the tail behavior of the distribution function depends on the number to of polled agents. 

Figure [IT] shows the behavior of the autocorrelation function for to = 60 and to = 100, with 
the same values of the other parameters phb = Pbh = 0.72 and p^h = Pbb = 0.85. For to = 100, the 
presence of the weakly unstable orbits is felt much stronger, which is reflected in 1) a very strong periodic 
component of the correlation function and 2) its slow decay. Even for to = 60, the correlation function 
does not decay fast enough compared to the typical duration of speculative bubbles to be in quantitative 
agreement with empirical data. This anomalously large correlation of the returns is obviously related to 
the deterministic dynamics of the returns. We thus expect that including stochastic noise due to a finite 
number N of agents (see below) and adding external noise due to "news" will whiten r t significantly. 

Figure [j~2| compares the correlation function for the returns time series r t oc p t — 1/2 and the volatility 
time series defined as |r t |. The volatility is an important measure of risks and thus plays an important 
role in portfolio managements and option pricing and hedging. Note that taking the absolute value of the 
return removes the one source of irregularity stemming from the change of sign of r t oc p t — 1/2 to focus 
on the local amplitudes. We observe in figure |l2| a significantly longer correlation time for the volatility. 
Moreover, the correlation function of the volatility first decays exponentially and then as a power law. 
This behavior has previously been documented in many econometric works (Ding et al. (1993), Ding 
and Granger (1996), Miiller et al. (1997), Dacorogna et al. (1998), Arneodo et al. (1998), Ballocchi et 
al. (1999), Muzy et al. (2001)). 
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6 Asymmetric cases 



We have seen that the symmetric case phb — Pbh an d phh — Pbb is plagued by the weakly unstable 
periodic orbits which put a strong and unrealistic imprint on the statistical properties of the return time 
series. It is natural to argue that breaking the symmetry will destroy the strength of these periodic orbits. 

From a behavioral point of view, it is also quite clear that the attitude of an investor is not symmetric. 
One can expect a priori a stronger bullish buy-and-hold index p^ — 1/2 than bearish buy-and-hold 
index pbh — 1/2: one is a priori more prone to hold a position in a bullish market than in a bearish 
one. Similarly, we expect a smaller bullish reversal index 1 — p^ than bearish reversal index 1 — p^: 
speculative bubbles are rarely seen on downward trends as it is much more common that increasing 
prices are favorably perceived and can be sustained much longer without reference to the fundamental 
price. 

Such an asymmetry has been clearly demonstrated empirically in the difference between the rate of 
occurrence and size of extreme drawdowns compared to drawups in stock market time series (Johansen 
and Sornette (2001)). Drawdowns (drawups) are defined as the cumulative losses (gains) from the last 
local maximum (minimum) to the next local minimum (maximum). Drawdowns and drawups are very 
interesting because they offer a more natural measure of real market risks than the variance, the value- 
at-risk or other measures based on fixed time scale distributions of returns. For the major stock market 
indices, there are very large drawdowns which are "outliers" while drawups do not exhibit such drastic 
change of regime. For major companies, drawups of amplitude larger than 15% occur at a rate about 
twice as large as the rate of drawdowns, but with lower absolute amplitude. 

Figure 13 compares the dynamics for the symmetric system (upper panel (a)) and for the asymmetric 
system (lower panel (b)). It is clear that, as expected, the number of periodic orbits decreases significantly 
in the asymmetric system. However, there are still an unrealistic number of negative bubbles. It is not 
possible to increase the asymmetry sufficiently strongly without exiting from the chaotic regime. This 
unrealistic feature is thus an intrinsic property and limitation of the present model. We shall indicate in 
the conclusion possible extensions and remedies. 

Figure [l4| compares the cumulative distributions of p — 1/2 for m = 60 for the symmetric and 
asymmetric cases. The strong effect of the weakly unstable periodic orbits observed in the periodic case 
has disappeared. In addition, the tail of the distribution decays faster in the asymmetric case, in better 
(but still not very good) agreement with empirical data. 

Figure |l5| shows the correlation function of the returns for a symmetric and an asymmetric case. In 
the asymmetric case, there is no trace of oscillations but the decay is slightly slower. 



7 Finite size effects 

Until now, our analysis has focused on the limit of an infinite number N — > oo of agents, in which 
each agent polls randomly m agents among N. In this limit, we have shown that, for a large domain 
in the parameter space, the dynamics of the returns is chaotic with interesting and qualitatively realistic 
properties. 

7.1 Finite-size effects in other models 

We now investigate finite-size effects resulting from a finite number N of interacting agents trading on 
the stock market. This issue of the role of the number of agents has recently been investigated vigor- 
ously with surprising results. First, Egenter et ai, (1999) studied the iV-dependence of the dynamical 
properties of price time series of the Kim-Markowitz (1989) and of the Lux-Marchesi (1999) models. 
They found that, if this number N goes to infinity, nearly periodic oscillations occur and the statistical 
properties of the price time series become completely unrealistic. Stauffer (1999) reviewed this work 
and others such as the Levy-Levy-Solomon (1995, 2000) model: realistically looking price fluctuations 
are obtained for N ~ 10 2 , but for N ~ 10 6 the prices vary smoothly in a nearly periodic and thus unre- 
alistic way. The model proposed by Farmer (1998) suffers from the same problem: with a few hundred 
investors, most investors are fundamentalists during calm times, but bursts of high volatility coincide 
with large fractions of noise traders. When N becomes much larger, the fraction of noise traders goes to 
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zero in contradiction to reality. On a somewhat different issue, Huang and Solomon (2001) have stud- 
ied finite-size effects in dynamical systems of price evolution with multiplicative noise. They find that 
the exponent of the Pareto law obtained in stochastic multiplicative market models is crucially affected 
by a finite ./V and may cause in the absence of an appropriate social policy extreme wealth inequality 
and market instability. Another model (apart from ours) where the market may stay realistic even for 
N — > oo seems to be the Cont-Bouchaud percolation model (2001). However, this only occurs for an 
unrealistic tuning of the percolation concentration to its critical value. Thus, in most cases, the limit 
N — ► oo leads to a behavior of the simulated markets which becomes quite smooth or periodic and thus 
predictable, in contrast to real markets. Our model which remains (deterministically) chaotic is thus a 
significant improvement upon this behavior. We trace this improvement on the highly nonlinear behav- 
ior resulting from the interplay between the imitative and contrarian behavior. It has thus been argued 
(Stauffer (1999)) that, if these previous models are good descriptions of markets, then real markets with 
their strong random fluctuations are dominated by a rather limited number of large players: this amounts 
to assume that the hundred most important investors or investment companies have much more influence 
than the millions of less wealthy private investors. 

There is another class of models, the minority games (Challet and Zhang (1997)), in which the dy- 
namics remains complex even in the limit N — > oo. It has been established that the fluctuations of 
the sum of the aggregate demand have an approximate scaling with similar sized fluctuations (volatil- 
ity/standard deviation) for any N and m for the scale scaled variable 2 m /N, where m is the memory 
length (Challet et al. (2000)). In a generalization, the so-called Grand Canonical version of the Mi- 
nority Game (Jefferies et al. (2001)), where the agents have a confidence threshold that prevents them 
from playing if their strategies have not been successful over the last T turns, the dynamics can depend 
more sensitively on N: as ./V becomes small, the dynamics can become quite different. For large N, the 
complexity remains. 

The difference between the limit N — ► oo considered up to now in this paper and the case of finite N 
is that p t is no more the fraction of bullish agents. For finite N, pt must be interpreted as the probability 
for an agent to be bullish. Of course, in the limit of large N, the law of large numbers ensures that the 
fraction of bullish agents becomes equal to the probability for an agent to be bullish. There are several 
ways to implement a finite-size effect. We here discuss only the two simplest ones. 



7.2 Finite external sampling of an infinite system 

Consider a system with an infinite number of agents for which the fraction p t of bullish agents is governed 
by the deterministic dynamics (^). At each time step t, let us sample a finite number N of them to 
determine the fraction of bullish agents. We get a number n, which is in general close but not exactly 
equal to Np t due to statistical fluctuations. The probability to find n bullish agents among N agents is 
indeed given by the binomial law 

Pr(n)=f^p n (l-p) JV - n . (25) 

This shows that the observed proportion fi = n/N of bullish agents is asymptotically normal with mean p 
and standard deviation — p)N : Pr(p) ~ Af(p, 1/ y/p{l — p)N). Iterating the sampling among 

N agents at each time step gives a noisy dynamics fit shadowing the true deterministic one. 

Figure |l6] compares the dynamics of the deterministic pt corresponding to N — > oo (panel (a)) with 
fit for a number N — m + 1 = 61 of sampled agents among the infinite ensemble of them (panel (b)). 
Panel (c) is the "noise" time series defined as p t —pt, i.e, by subtracting the time series of panel (a) from 
the time series of panel (b). The noise time series of panel (c) thus represents the statistical fluctuations 
due to the finite sampling of agents 'opinions. Figure |l6|-b shows the characteristic volatility clusters 
which is one of the most important stylized properties of empirical time series. 

For large N, we can write 

fit = p t + 1 W t (26) 

y/p t (l-pt)N 

where {Wt} are iid gaussian variables with zero mean and unit variance. Therefore, the correlation 
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function Corr n(t) at lag r ^ is obtained from that for N — > oo by multiplication by a constant factor: 

Cott n (t) = rr^TI an/ ? ^ xCorr °o(r) and r ^ 0, (27) 

E[l/{p(l-p)}] +A^Var(p) 

~ Corr^r) for larger, (28) 

where E[x] denotes the expectation of a; with respect to the continuous invariant measure of the dynami- 
cal system (0). Note that E[l / {p(l — p)}] always exists for m < oo since the support of the continuous 
measure of (g) with respect to Lebesgue measure is bounded from below by a value strictly larger than 
and from above by a value strictly less than 1 . Figure [n| shows that the correlation function of p t is very 
close to that of the deterministic trajectory p t . 

To quantify further the impact of the statistical noise stemming from the finite size of the market, 



figures [18j and |19| show the return maps of p t , i.e, pt+i as a function of pt, for to = 60 polled agents 
among a total number N = 61 of agents (fig H)) and N = 600 (fig |9|). 

Figure ^ shows the price trajectory obtained by 7r t = 7r t _i exp[f t ] in linear and logarithmic scale. 
The super-exponential acceleration of the price giving rise to sharp peaks in the semi-logarithmic repre- 
sentation (Roehner and Sornette (1998)) is clearly visible. 



7.3 Finite number N of agents 

We now introduce a genuine finite stock market with N agents. We assume that the agents do not know 
the exact number N of agents in the market (this is realistic) and they are in contact with only to other 
agents that they poll at each time period. Not knowing the true value of N but assuming it to be large, 
it is rational for them to develop the best predictor of the dynamics by assuming the ideal case of an 
infinite number of agents with to polled agents and thus use the deterministic dynamics (Q) as their best 
predictor. 

At each time period t, each agent thus chooses randomly m agents that she polls. She then counts 
the number of bullish and bearish agents among her polled sample of m agents. This number divided by 
to gives her an estimation p t of the probability p t be to bullish at time t. Introducing this estimation in 
the deterministic equation (||), the agent obtains a forecast p' of the true probability p' to be bullish at the 
next time step. 

Results of the simulations of this model are shown in figure [HJ We observe a significantly stronger 
"noise" compared to the previous section, which is expected since the noise is itself injected in the 
dynamical equation at each time step. As a consequence, the correlation function of the returns and 
of the volatility decay faster than their deterministic counterpart. The correlation of the volatility still 
decays about ten times slower than the correlation of the returns, but this clustering of volatility is not 
sufficiently strong compared to empirical facts. 

Other more realistic models of a finite number of agents can be introduced. For instance, at time t, 
consider an agent among the N. She chooses to other agents randomly and polls them. Each of them is 
either bullish or bearish as a result of decisions taken during the previous time period. She then counts the 
number of bullish agents among the to, and then determines her new attitude using the rules ([[]). If she is 
polled at time t + 1 by another agent, her attitude will be the one determined from t to t+ 1. In this way, 
we never refer to the deterministic dynamics p t but only to its underlying rules. As a consequence, this 
deterministic dynamics does not exert an attraction that minimizes the effect of statistical fluctuations 
due to finite sizes. This approach is similar to going from a Fokker-Planck equation (equation (^)) to a 
Langevin equation with finite-size effects. This class of models will be investigated elsewhere. 



8 Conclusions 

The traditional concept of stock market dynamics envisions a stream of stochastic "news" that may move 
prices in random directions. This paper, in contrast, demonstrates that certain types of deterministic 
behavior — mimicry and contradictory behavior alone — can already lead to chaotic prices. 

If the market prices are assumed to follow the p t behavior, our description refers to the well-known 
evolution of the speculative bubbles. Such apparent regularities often occur in the stock market and form 
the basis of the so-called "technical analysis" whereby traders attempt to predict future price movements 
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by extrapolating certain patterns from recent historical prices. Our model provides an explanation of 
birth, life and death of the speculative bubbles in this context. 

While the traditional theory of rational anticipations exhibits and emphasizes self -reinforcing mech- 
anisms, without either predicting their inception nor their collapse, the strength of our model is to justify 
the occurrence of speculative bubbles. It allows for their collapse by taking into account the combination 
of mimetic and antagonistic behavior in the formation of expectations about prices. 

The specific feature of the model is to combine these two Keynesian aspects of speculation and 
enterprise and to derive from them behavioral rules based on collective opinion: the agents can adopt an 
imitative and gregarious behavior, or, on the contrary, anticipate a reversal of tendency, thereby detaching 
themselves from the current trend. It is this duality, the continuous coexistence of these two elements, 
which is at the origin of the properties of our model: chaotic behavior and the generation of bubbles. 

It is a common wisdom that deterministic chaos leads to fundamental limits of predictability because 
the tiny inevitable fluctuations in those chaotic systems quickly snowball in unpredictable ways. This 
has been investigated in relation with for instance long-term weather patterns. However, in the context of 
our models, we have shown that the chaotic dynamics of the returns alone cannot be the limiting factor 
for predictability, as it contains too much residual correlations. Endogenous fluctuations due to finite- 
size effects and external news (noise) seem to be needed as important factors leading to the observed 
randomness of stock market prices. The relation between these extrinsic factors and the intrinsic ones 
studied in this paper will be explored elsewhere. 

Remark and Acknowledgements: This paper is an outgrowth and extension of unpublished work by 
three of us (AC, JPE, AM) which was in turn based on the Ph.D. of Anne Corcos. We are grateful to J.V. 
Andersen for useful discussions. This work was partially supported by the Fonds National Suisse (JPE 
and AM) and by the James S. Mc Donnell Foundation 21st century scientist award/studying complex 
system (DS). 
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Appendix 

We expand F m (p) around the fixed point p = 1/2, so that, using the symmetry of F m {p) 

F m (p) = \ + F' m {\/2) -(p-1)+ ^ 3) (V2) • (p £f + • • • 

First of all, it is obvious to show by recursion that 

F m (l/2) = l-2 5m (l/2)- g' m (l/2) 
F ^+i) {1/2) = _ 2{2k+ 1)^(1/2) - 5 (f +i )(1/2) if fc > 

The problem thus amounts to calculating the derivatives of g m . 
Some simple algebraic manipulations allow to obtain 



0m (?) = m E 

m — 1 

= - m E 

J=0 



TO — 1 
j 

m — 1 



p rol - j (i-p) j [/(^J-/ 

p ml -^(l-f^A 1 / m (j), 



i + i 



TO 



where Ai/ TO (-) is the first order discrete derivative of / ( — ), which yields 



By recursion, it is easy to prove that 

7 (fe) 



m— 1 

tE 

3=0 



TO — 1 
j 



Ai/ m (j) . 



/ , ^ ■ m — fc / i 



2 m-fc fc! 



j=0 



and Afe/ m (-) is the fc" 1 order discrete derivative of / (— ): 



A fe / ro (i) = ^( )(-l)7 



i=0 



j + i 



Finally, 



2 m-2fe-l ( 2 fc)! 



-(2fc + l) ( m . 2fc )A 2fe / ro (j) 
j=0 \ j / 



^ m— 2fc — 1 

2fc + l S 



to - 2fc - 1 



A2fc+l/r? 
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Figure 3: The time series for the same parameter values as in Fig. 2. Note that, for phh = Pbb = 0.75, one 
has convergence to a bullish equilibrium, for 0.76 a bullish period 2, for 0.77 a bullish, but chaotic behavior. 
The most interesting case is p^h = Pbb = 0.85, where calm periods alternate in a seemingly random fashion 
with speculative bubbles. 
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Figure 4: Evolution of the system over 10000 time steps for N = oo, m = 60 polled agents and the 
parameters p hb = p bh = 0.72 and p hh = p bb = 0.85. 
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Figure 6: The logarithm of F m (p) — p versus the logarithm of p — 1/2 for three different values of 
m = 30, 60 and 100, with pht = Pbh = 0.72 and p^h = Pbb = 0.85. Note the transition from a slope 1 to a 
large effective slope before the reinjection due to the contrarian mechanism. 
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Figure 7: Approximation of the function F m {p) — \ by the function f(p) = [1 + a] (p + |j + (3 [p + |j 
over different p-intervals, for m = 60 interacting agents and parameters p^ = p^ = 0.72 and p^h = Pbb = 
0.85. 
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Figure 8: ( pt _\/2)' 2 versus * to qualify the finite time singularity predicted by ( |23| ) for m = 60 with 
Phb = Pbh = 0.72 and p^ = = 0.85. The points are obtained from the time series pt and the straight 
continuous lines are the best linear fits. The horizontal dashed lines indicate the levels p — 1/2 = 0.05, 0.01 
and 0.2 to demonstrate that most of the visited values are close to the unstable fixed point, which determines 
the effective value of the nonlinear exponent p 3. 
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Figure 9: Cumulative distribution for m = 60 polled agents and the parameters phb = Pbh = 0.72 and 
Phh = Pbb = 0.85. 
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Figure 11: The upper panels represent the time series p t for m = 60 (left) and m = 100 (right). The 
lower panels represents the corresponding autocorrelation function of r t oc p — 1/2 for m = 60 (left) and 
m = 100 (right) with the same parameters = = 0.72 and phh = Pbb = 0.85. 
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Figure 13: Time evolution of pt over 10000 time steps for m = 60 polled agents in (a) a symmetric case 
Phb = Pbh = 0.72 and p hh = p bb = 0.85 and (b) an asymmetric case p hb = 0.72, p bh = 0.74, p hh = 0.85 
and p bb = 0.87. 
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Figure 14: Distribution function of p — 1/2 for m = 60 polled agents and parameters phb = Pbh = 0.72 
and phh = Pbb = 0.85 (dashed line) and p^b = 0.72, pbh = 0.74, phh = 0.85 and pbb = 0.87 (continuous 
line). 
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ure 15: Correlation function for m = 60 polled agents and parameters p^b = Pbh = 0.72 and p^h 
= 0.85 (dashed line) and phb = 0.72, pbh = 0.74, phh = 0.85 and pbb = 0.87 (continuous line). 
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Figure 16: Time evolution of pt over 10000 time steps for m = 60 polled agents with (a) N = oo, (b) 
N = m + 1 = 61 agents and parameters p^ = = 0.72 and phh = p^ = 0.85. The panel (c) represents 
the noise due to the finite size of the system and is obtained by subtracting the time series in panel (a) from 
the time series in panel (b). 
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Figure 17: Correlation function for m = 60 polled agents with N = oo (thin line), N = 600 (dashed line) 
and N = 61 (continuous line) agents and parameters phb = Pbh = 0.72 and phh = Pbb = 0.85. 
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Figure 18: Return map of the fraction of bullish agents for m = 60 polled agents among N = 61 agents 
(points) and the deterministic trajectory (continuous line) corresponding to N = oo agents. The parameters 
are p hb = p bh = 0.72 and p hh = p bb = 0.85. 
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Figure 19: Return map of the fraction of bullish agents for m = 60 polled agents among N = 600 agents 
(points) and the deterministic trajectory (continuous line) corresponding to N = oo agents. The parameters 
are p hb = p bh = 0.72 and p hh = p bb = 0.85. 
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Figure 20: Upper panel: return trajectory f% = jpt — 1/2 for m = 100, N = 100, phb = Pbh = 0.72 and 
Phh = Pbb = 0.85 and 7 = 0.01. Middle panel: price trajectory obtained by ir t = wt-l exp[f t ] . Lower 
panel: same as the middle panel with it t shown in logarithmic scale. Note the "flat trough-sharp peak" 
structure of the log-price trajectory (Roehner Sornette (1998)). 
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Figure 21: Evolution of the system over 10000 time steps for m = 60 polled agents with (upper panel) 
N = oo, (second panel) N = m + 1 = 61 and parameters p^ = p^ = 0.72 and p^h = Pbb = 0-85. The 
lower panel represents the "noise" introduced by the finite size of the system and is obtained by subtracting 
the upper panel from the second panel. 
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